Transcriptomic Analysis of Papillary Thyroid Cancer: A Focus on Immune-Subtyping, Oncogenic Fusion, and Recurrence

Objectives@#. Thyroid cancer is the most common endocrine tumor, with rapidly increasing incidence worldwide. However, its transcriptomic characteristics associated with immunological signatures, driver fusions, and recurrence markers remain unclear. We aimed to investigate the transcriptomic characteristics of advanced papillary thyroid cancer. @*Methods@#. This study included 282 papillary thyroid cancer tumor samples and 155 normal samples from Chungnam National University Hospital and Seoul National University Hospital. Transcriptomic quantification was determined by high-throughput RNA sequencing. We investigated the associations of clinical parameters and molecular signatures using RNA sequencing. We validated predictive biomarkers using the Cancer Genome Atlas database. @*Results@#. Through a comparison of differentially expressed genes, gene sets, and pathways in papillary thyroid cancer compared to normal tumor-adjacent tissue, we found increased immune signaling associated with cytokines or T cells and decreased thyroid hormone synthetic pathways. In addition, patients with recurrence presented increased CD8+ T-cell and Th1-cell signatures. Interestingly, we found differentially overexpressed genes related to immune-escape signaling such as CTLA4, IDO1, LAG3, and PDCD1 in advanced papillary thyroid cancer with a low thyroid differentiation score. Fusion analysis showed that the PI3K and mitogen-activated protein kinase (MAPK) signaling pathways were regulated differently according to the RET fusion partner genes (CCDC6 or NCOA4). Finally, we identified HOXD9 as a novel molecular biomarker that predicts the recurrence of thyroid cancer in addition to known risk factors (tumor size, lymph node metastasis, and extrathyroidal extension). @*Conclusion@#. We identified a high association with immune-escape signaling in the immune-hot group with aggressive clinical characteristics among Korean thyroid cancer patients. Moreover, RET fusion differentially regulated PI3K and MAPK signaling depending on the partner gene of RET, and HOXD9 was found to be a recurrence marker for advanced papillary thyroid cancer.

Bioinformatic identification of prognostic signature defined by copy number alteration and expression of CCNE1 in non-muscle invasive bladder cancer

Non-muscle invasive bladder cancer (NMIBC) patients frequently fail to respond to treatment and experience disease progression because of their clinical and biological diversity. In this study, we identify a prognostic molecular signature for predicting the heterogeneity of NMIBC by using an integrative analysis of copy number and gene expression data. We analyzed the copy number and gene expression profiles of 404 patients with bladder cancer obtained from The Cancer Genome Atlas (TCGA) consortium. Of the 14 molecules with significant copy number alterations that were previously reported, 13 were significantly correlated with copy number and expression changes. Prognostic gene sets based on the 13 genes were developed, and their prognostic values were verified in three independent patient cohorts (n=501). Among them, a signature of CCNE1 and its coexpressed genes was significantly associated with disease progression and validated in the independent cohorts. The CCNE1 signature was an independent risk factor based on the result of a multivariate analysis (hazard ratio=6.849, 95% confidence interval=1.613–29.092, P=0.009). Finally, gene network and upstream regulator analyses revealed that NMIBC progression is potentially mediated by CCND1-CCNE1-SP1 pathways. The prognostic molecular signature defined by copy number and expression changes of CCNE1 suggests a novel diagnostic tool for predicting the likelihood of NMIBC progression.

A polymorphic minisatellite region of BORIS regulates gene expression and its rare variants correlate with lung cancer susceptibility

Aberrant expression of BORIS/CTCFL (Brother of the Regulator of Imprinted Sites/CTCF-like protein) is reported in different malignancies. In this study, we characterized the entire promoter region of BORIS/CTCFL, including the CpG islands, to assess the relationship between BORIS expression and lung cancer. To simplify the construction of luciferase reporter cassettes with various-sized portions of the upstream region, genomic copies of BORIS were isolated using TAR cloning technology. We analyzed three promoter blocks: the GATA/CCAAT box, the CpG islands and the minisatellite region BORIS-MS2. Polymorphic minisatellite sequences were isolated from genomic DNA prepared from the blood of controls and cases. Of the three promoter blocks, the GATA/CCAAT box was determined to be a critical element of the core promoter, while the CpG islands and the BORIS-MS2 minisatellite region were found to act as regulators. Interestingly, the polymorphic minisatellite region BORIS-MS2 was identified as a negative regulator that repressed the expression levels of luciferase reporter cassettes less effectively in cancer cells compared with normal cells. We also examined the association between the size of BORIS-MS2 and lung cancer in a case–control study with 590 controls and 206 lung cancer cases. Rare alleles of BORIS-MS2 were associated with a statistically significantly increased risk of lung cancer (odds ratio, 2.04; 95% confidence interval, 1.02–4.08; and P=0.039). To conclude, our data provide information on the organization of the BORIS promoter region and gene regulation in normal and cancer cells. In addition, we propose that specific alleles of the BORIS-MS2 region could be used to identify the risk for lung cancer.

A Database of Gene Expression Profiles of Korean Cancer Genome

Because there are clear molecular differences entailing different treatment effectiveness between Korean and non-Korean cancer patients, identifying distinct molecular characteristics of Korean cancers is profoundly important. Here, we report a web-based data repository, namely Korean Cancer Genome Database (KCGD), for searching gene signatures associated with Korean cancer patients. Currently, a total of 1,403 cancer genomics data were collected, processed and stored in our repository, an ever-growing database. We incorporated most widely used statistical survival analysis methods including the Cox proportional hazard model, log-rank test and Kaplan-Meier plot to provide instant significance estimation for searched molecules. As an initial repository with the aim of Korean-specific marker detection, KCGD would be a promising web application for users without bioinformatics expertise to identify significant factors associated with cancer in Korean.

Pro-oncogenic potential of NM23-H2 in hepatocellular carcinoma

NM23 is a family of structurally and functionally conserved proteins known as nucleoside diphosphate kinases (NDPK). There is abundant mRNA expression of NM23-H1, NM23-H2, or a read through transcript (NM23-LV) in the primary sites of hepatocellular carcinoma (HCC). Although the NM23-H1 protein is implicated as a metastasis suppressor, the role of NM23-H2 appears to be less understood. Thus, the aim of this study was to examine whether NM23-H2 is associated with hepatocarcinogenesis. The level of NM23-H2 expression in tumor tissues and the surrounding matrix appeared to be independent of etiology and tumor differentiation. Its subcellular localization was confined to mainly the cytoplasm and to a lesser extent in the nucleus. Ectopic expression of NM23-H2 in NIH3T3 fibroblasts and HLK3 hepatocytes showed a transformed morphology, enhanced focus formation, and allowed anchorage-independent growth. Finally, NIH3T3 fibroblasts and HLK3 hepatocytes stably expressing NM23-H2 produced tumors in athymic mice and showed c-Myc over-expression. In addition, NF-kappaB and cyclin D1 expression were also increased by NM23-H2. Lentiviral delivery of NM23-H2 shRNA inhibited tumor growth of xenotransplanted tumors produced from HLK3 cells stably expressing NM23-H2. Collectively, these results indicate that NM23-H2 may be pro-oncogenic in hepatocarcinogenesis.

Genes Associated with Recurrence of Hepatocellular Carcinoma: Integrated Analysis by Gene Expression and Methylation Profiling

Gene expression is suppressed by DNA methylation. The goal of this study was to identify genes whose CpG site methylation and mRNA expression are associated with recurrence after surgical resection for hepatocellular carcinoma (HCC). Sixty-two HCCs were examined by both whole genome DNA methylation and transcriptome analysis. The Cox model was used to select genes associated with recurrence. A validation was performed in an independent cohort of 66 HCC patients. Among fifty-nine common genes, increased CpG site methylation and decreased mRNA expression were associated with recurrence for 12 genes (Group A), whereas decreased CpG site methylation and increased mRNA expression were associated with recurrence for 25 genes (Group B). The remaining 22 genes were defined as Group C. Complement factor H (CFH) and myosin VIIA and Rab interacting protein (MYRIP) in Group A; proline/serine-rich coiled-coil 1 (PSRC1), meiotic recombination 11 homolog A (MRE11A), and myosin IE (MYO1E) in Group B; and autophagy-related protein LC3 A (MAP1LC3A), and NADH dehydrogenase 1 alpha subcomplex assembly factor 1 (NDUFAF1) in Group C were validated. In conclusion, potential tumor suppressor (CFH, MYRIP) and oncogenes (PSRC1, MRE11A, MYO1E) in HCC are reported. The regulation of individual genes by methylation in hepatocarcinogenesis needs to be validated.

Bioinformatics Resources of the Korean Bioinformation Center (KOBIC)

The Korean Bioinformation Center (KOBIC) is a national bioinformatics research center in Korea. We developed many bioinformatics algorithms and applications to facilitate the biological interpretation of OMICS data. Here we present an introduction to major bioinformatics resources of databases and tools developed at KOBIC. These resources are classified into three main fields: genome, proteome, and literature. In the genomic resources, we constructed several pipelines for next generation sequencing (NGS) data processing and developed analysis algorithms and web-based database servers including miRGator, ESTpass, and CleanEST. We also built integrated databases and servers for microarray expression data such as MDCDP. As for the proteome data, VnD database, WDAC, Localizome, and CHARMM_HM web servers are available for various purposes. We constructed IntoPub server and Patome database in the literature field. We continue constructing and maintaining the bioinformatics infrastructure and developing algorithms.

Identification of molecular markers for the oncogenic differentiation of hepatocellular carcinoma

The aim of this study was to identify molecular markers associated with oncogenic differentiation in hepatocellular carcinoma (HCC). Using an unsupervised clustering method with a cDNA microarray, HCC (T) gene expression profiles and corresponding non-tumor tissues (NT) from 40 patients were analyzed. Of total 217 genes, 72 were expressed preferentially in HCC tissues. Among 186 differentially regulated genes, there were molecular chaperone and tumor suppressor gene clusters in the Edmondson grades I and II (GI/II) subclass compared with the liver cirrhosis (LC) subclass. The Edmondson grades III and IV (GIII/IV) subclass with a poor survival (P = 0.0133) contained 122 differentially regulated genes with a cluster containing various metastasis- and invasion-related genes compared with the GI/II subclass. Immunohistochemical analysis revealed that ANXA2, one of the 72 genes preferentially expressed in HCC, was over-expressed in the sinusoidal endothelium and in malignant hepatocytes in HCC. The genes identified in the HCC subclasses will be useful molecular markers for the genesis and progression of HCC. In addition, ANXA2 might be a novel marker for tumor angiogenesis in HCC.